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Abstract 

Solution of the nonlinear Klein-Gordon equation perturbed by small 
external force is investigated. The frequency of perturbation varies slowly 
and passes through a resonance. The resonance generates a solitary pack- 
ets of waves. Full asymptotic description of this process is presented. 

Introduction 

This work is devoted to the problem on a generation of solitary packets 
of waves by a small external driving force. We propose a new approach for 
generation of solitary packets of waves. We demonstrate that for perturbed 
nonlinear Klein-Gordon equation. 

In our approach the wave packets appear due to passing of external driving 
force through resonance. After the resonance the envelope function of the wave 
packet is determined by nonlinear Schrodinger equation (NLSE). In the most 
important cases the envelope function is a sequence of solitary waves which are 
called solitons. The wave packets with the solitons as the envelope function are 
propagated without a deformation. The parameters of the solitons are obviously 
defined by the value of the driving force on a resonance curve. 

Here we give the mathematical basis for the proposed approach. This basis 
allows to derive explicit formulas which define parameters for the solitary packets 
of waves with respect to the external driving force. Generation of the solitary 
packets of waves by a small driving force is described in detail. The formulas for 
the asymptotic solution before, after and in the neighborhood of the resonance 
curve are obtained. 

Proposed approach is based on a local resonance phenomenon. The local 
resonance in linear ordinary differential equations was investigated in papers 
|21 E]. Later this phenomenon was investigated in partial differential equations 
in linear case 0] and in weak nonlinear case |SJ|S]. As was shown in these papers 
the amplitude of the wave which crosses the local resonance increases by linear 
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way. The increase of the amplitude is proportional to the width of the local 
resonance layer. 

After resonance a special proportion between the order of solution and scales 
of independent variables appears. This magic proportion gives NLSE for enve- 
lope function. The deriving of NLSE in such case is well known [3 |SJ |S] and 



Important kind of solution of NLSE is solitons. It is known the phenomenon 
of solitary waves generation for some nonlinear equations due to the modulation 
instability For example, there exist the detailed analytical description of 
disappearance (generation) of soliton due to the modulation instability in the 
case of Kadomtsev-Petviashvili equation ^21- Some results on soliton appear- 
ance in nonlinear Schrodinger equation due to instability are presented in |13|. 
It is known the structural instability of the solitons for Davey-Stewartson equa- 
tion ^3]. Perturbations used the different kinds of instability of the solution do 
not allow to obtain the solitons with given parameters. 

The resonance generation of solitary waves by small external force is known 
due to computer simulation 1 . Although computer simulations justify the 
possibility of soliton generation by an external driving force but also do not 
allow to connect the the parameters of the soliton and perturbation. Therefore 
the problem on soliton generation with the given parameters was still open. 

The goal of this paper is the following: to demonstrate that the process 
of solitary waves generation due to local resonance is universal. This process 
allows to control parameters of generated waves. Earlier this phenomenon in 
the case of nonlinear Schrodinger equation was asymptotically investigated in 
|15| . In this work we consider the similar phenomenon in the nonlinear Klein- 
Gordon equation. Our approach demonstrates that solitary waves with the 
given parameters can be obtained for nonlinear wave systems. 

This paper has the following structure. The first section contains the main 
result and example. The second section contains the asymptotic construction 
in the pre-resonance domain. In the third section we construct the asymptotic 
solution in the neighborhood of the resonance curve. The fourth section of 
the paper is devoted to construction of the post-resonance asymptotics. All 
asymptotics are matched. 

1 Main result 

Let us consider the Klein-Gordon equation with cubic nonlinearity 



d 2 U - d 2 x U + U + jU 3 = £ 2 f(£x)exp\i {£ y X > \ +c.c, 0<e<l. (1) 



Here 7 = const; f(y) is smooth and rapidly vanish as y — > ±00. The phase 
function S(y, z) of the driving force and all derivates with respect to y, z are 
bounded. 

We construct the formal asymptotic solution of the WKB-type using the 
combination of method of multiple scales ^B] and matching method ^Zj . Below 



justified HJ3]. 
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we use scaled variables 

j ' — ' & } ^3 — t ^ — 1 j 2 ■ 

The leading-order term of the asymptotic solution has an order e 2 and os- 
cillates with frequency of the driving force. The resonance curve is determined 

by 

i(t 2 ,x 2 ) = (d t2 s) 2 -(d X2 s) 2 -i = o. 

We assume that the curve S = const does not touch the resonance curve: 

d X2 w X2 s - d t2 w t2 s ^ o. 

The frequency of the forced oscillations and frequency of the eigen oscillations 
of the linearized Klein- Gordon equation are equal on the resonance curve. It 
yields the local resonance layer in the neighborhood of the curve l(x 2 ,t 2 ) = 0. 
The asymptotics of the WKB-type is not valid into this layer. The leading- 
order term of the asymptotics is defined by Fresnel integral. After resonance 
the solution has an order e and oscillates. The envelope function satisfies NLSE. 
The accurate formulation of the result for this paper is following 

Theorem 1 In the domain — I 3> s the asymptotic solution of 0) modulo 
0(e N+1 ) has the form 

N 

U = Y,£ n U(t,x,s), 

where the leading- order term is 
2 f 

U= - j exp(iS{t 2 , x 2 )/e 2 ) + c.c. 

The higher-order terms are determined from the recurrent system of algebraic 
equations |3J), Q). 

In the domain \l\ <§C 1 the asymptotic solution of modulo 0(e N+1 ) has the 
form 

N 

En 
e n W (t u x 1 ,t 2 ,x 2 ,e). 

n>l 

The leading-order term is 

1 1 

W (ti,xi,t 2 ,x 2 ,e) =Wi exp(iS(t 2 ,x 2 )/e ) + c.c. 
1 

The function Wi is solution of the equation 

2id t2 Sd tl Wi -2id X2 Sd xl Wi -A Wi= f 

with the zero condition as A — > —00, where X(xi,ti,e) — l(t 2 ,x 2 )/e. The higher- 
order terms are either solutions of the problems for equations JJ7| ) with zero 
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conditions as A — > — oo or solutions of the algebraic equations flUl) . 

In the domain I e the asymptotic solution of £7J) modulo 0(e N+1 ) has the 

form 

N n-2 / , 



U(x,t,e) = J2 £n J2 lnk ^( ^exp{±i$(a; 2 ,t 2 )/£ 2 } $ ±*(a:i,ti,t2 

n=l fc=0 ^ ±* 

+ 2^ exp{ix(a;2,i2)/e } * x (xi,ti,t 2 ) 

T/ie phase function $ satisfies eikonal equation 

(d t2 $) 2 -(d X2 $) 2 -l=0 

wii/i conditions 

$|, =0 =S|/= , d t2 $| /=0 = d t2 S| ;=0 . 

TTie envelope function of the leading-order term is a solution of the nonlinear 
Schrodinger equation 

1,0 „ 1,0 1,0 1,0 1,0 

2idt 2 <£d t2 W $ + <9| * * + i[d t2 $ - <%: 2 $] * * + 7l * *l * * = 0, 



where the £ is defined from 



1,0 

T/ie initial condition for \j> $ is 

1,0 Z" 00 / >CT 

* $|i=o = / dcr/(a;i)exp(z / d/*A(a;i,ti,e)), 

./-oo JO 

Integration in this integral is realized in the line of characteristic direction con- 

n,k 

nected with \2ty . The coefficient \f $ is determined from Cauchy prob- 

n,k 

lems for linearized Schrodinger equation T7ie coefficients ^ x , X G fe 

are determined from algebraic equations fifyfl. The set K n ^ is defined by 

Ki t0 = ±&, K 2 , = ±<f>,±S, K 3 ,i = ±*, K ntk = Ua + f3 + 5, 

where 

aeK juh , 0€K hi i 2 , 5 € K h j 3 , ji+ h+ iz = n, h + h + h = k 
The set K' k = X„, fe \{±$}. 
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This theorem is a direct consequence of theorems 2,3 and 4 which are proved 
in the next parts of the paper. 

To illustrate the theorem we consider equation with the simplest of driv- 
ing force. Let the phase function of the driven be S = (s 2 t) 2 /2. In this case 
the curve of the local resonance is the line t% = 1. In the domain t 2 > 1 the 
leading-order term of the asymptotics satisfies the Cauchy problem 

1,0 1,0 1,0 1,0 

2id t2 * * + ^ * * + 7l * *l * * = 0, 
**|t a =i=/(0(l + i)V5F- 

The solution of this Cauchy problem contains solitary waves if the initial data 
is sufficiently large |12) . 

Remark on WKB asymptotics. Theoremndescribes the special asymp- 
totic solution of equation (JTJ. It is defined by the driving force. One can add 
any solution of the WKB-type ^B] of the order e 2 to this constructed solution. 
It leads to an asymptotic solution for equation of the form 

N 

U = U{t, x, e) + e" U [t, x, e). 

n>2 

n 

The coefficients U (t, x, e) of the asymptotics are calculated by standard meth- 
ods of WKB-theory. This additional term leads to ponderous formulas and does 
not change the leading-order term of the asymptotics constructed in theorem ^ 



2 Pre-resonance expansion 

In this section the formal asymptotic solution is constructed in the domain 
before resonance. The asymptotic expansion has the form of the WKB-type. 
The leading-order term of the asymptotics has the order of the driving force and 
oscillates with its frequency. The constructed asymptotics is valid as —I 3> s. 
The result of this section is formulated below. 

Let us construct the formal asymptotic solution for equation Q in the form 

N 

U = E £ " U e), (2) 

n>2 

where 

n ^-^ n f S(t 2 , X2) 1 

U= 2^ Uk {t 2 ,X2,£X)exjp i ik >. 

feen„ I e J 

Set f2 n for the higher-order term with the number n is described by the formula 




{±1}, n<5; 
{±l,±3,...,±(2/ + 3)}, Z=[(n-6)/4], n > 6. 
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The functions JJk and JJ-k are complex conjugated. 

Let us substitute ((2J) in equation and collect the terms of the same order 
of e. As a result we obtain a recurrent sequence of algebraic equations. 

2 f 

Ui=- L v (3) 
Ui=2i XlJ p X2 , (4) 
\2x2 



where 



4 2if[d t2 Sd t2 l - d X2 Sd X2 l] - A(d X2 S) 2 d 2 x J 

ui= Ta 

2id t Jd t2 S + d*J + idlSf 



i = (a t2 sy-(d X2 sy-i 



(5) 



The curve where the phase function S satisfies eikonal equation is called the 
resonance curve 

l[S] = (d t2 S) 2 - (0 X2 S) 2 - I = 0. (6) 

n 

The amplitude JJi has a singularity on this curve. 
The formula for the n-th order term has the form 



w l 



If n-4 n—2 n-2 n — 2 n-2 



d 2 2t2 Uk +2ikd t2 Sd t2 Uk +ikS t2t2 Uk -2ikd X2 Sd X2 Uk ~ikd 2 2 S U k 



n — 2 n — 3 n—4 n—1 ^ — ^ n± n 2 n 3 

d XlXl Uk-2d 2 XlX2 Uk-d 2 X2X2 Uk~2ikd X2 Sd Xl Uk+1 \^ U kl Uk 2 Uk 3 

£-'1+^2+^3 = ^ 

fceo„ 

(7) 

n 

Lemma 1 The coefficient JJk has the following behaviour 

Uk=o{r (n - k) ), fc>o, i^-o. (8) 

Proof. Let us prove this lemma as k = 1. The validity of formula © for 

n = 2,3,4 directly obtains from ©, Q), JSJ. Suppose now that this formula is 
n-l 

valid for the term JJ i . The increase of the order of the singularity as I — > 
takes place due to differentiation with respect to X2,t2 and the nonlinear term 
in formula Q . Differentiation of the terms in formula J7J) leads to formula (JSJ ■ 

n 

Let us consider JJk for k > 1. The validity of formula (JSJ) for small values 
of n and k obtains by direct calculations. Consider the n— th order term. It 
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contain the terms with different values of k. The higher-order term with k = 3 
have the greatest order of singularity. 

u s = o(r< n - 3 >), i^-o. (9) 

n-4 2 2 

It takes place because the right hand side of {7|) contains the term U ±i U±iU±i- 
The calculation of the order of singularity for this term leads to formula ©. 

n 3 ni iii 

The terms of the type of U ±3 U^k 1 U±k 1 , n\ + n 2 + m = n lead to weak 
singularities, for example for k\ = 3 we obtain the order of singularity is equal 
to 71 — 9. 

™! ™ 2 "3 . _ 

Consider nonlinear term UkiUk 2 Uk 3 from right hand side of J7J) when the 
number of the higher-order term is equal to n. Calculate the order of singularity 
for this term using the (n — 1)— th step of induction. Indexes of the amplitudes 
are connected by formulas 

n\ + n 2 + n 3 = n, k\ + fc 2 + k 3 = k. 

Using ijSJ for n±, n% < n we obtain that the order of the singularity for this 
term is equal to (n — k). 

The right hand side of (JJJ contains derivatives of previous terms with respect 
to X2,t 2 - It leads to increase of the order of the singularity but the leading order 
nevertheless we obtain from nonlinear terms. The lemma is proved. 

The domain of validity as I — > — for formal asymptotic solution in the form 
J2J is defined by 

e n+1 U 
=— < 1. 



e" U 

It yields 

-I > e. 

Using these lemmas we obtain the asymptotic representation for @ as i -> —0 

N 00 

U = £ " Yl exp{ikS/e 2 } I -»■ -0. (10) 

n=2 ken n j=-( n -fc) 

The following theorem is proved. 



Theorem 2 In the domain — / 3> e the formal asymptotic solution of equation 

n+l 

PP modulo 0(e N+1 ) has the form |5j) . The coefficients of the asymptotics JJ k 
are defined from algebraic equations ify) . ^j), (3|). TO. 



3 Internal asymptotics 

This part of the paper contains the asymptotic construction of the solution for 
equation Q in the neighborhood of the curve I — 0. The domain of validity 
of this asymptotics intersects with domain of validity of expansion © . These 
expansions are matched. 
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Theorem 3 In the domain \l\ -C 1 the formal asymptotic solution for equation 
modulo 0(e N+1 ) has the form 

N 

x — -\ n 

U = W (t 1 ,x 1 ,t 2 ,x 2 ,e), (11) 

n>l 

where 

n S. — > n I S(t 2 X2I I 

W= }^ Wk (x 2 ,t 2 ,x u ti)expUk 2 \, (12) 

The function Wk > k — \ is solution of the problem for equation with zero 
condition as X — > —00 and solutions of algebraic equations fiyp in the case 

n n 

k =/= 1. T/ie functions Wk and W—k are complex conjugated. 

There is an essential difference between asymptotics 111 II) and external pre- 
resonance asymptotics (J2J). First the leading order term in l|lljl has an order e 
in contrast the leading order term in (J2J has an order e 2 . Second the coefficients 
of asymptotics (f 1 If) depend on fast variables x\ = x 2 /e and t\ = t 2 je. 

The proof of theorem [3] consists in three steps. First we derive equations for 
coefficients of the asymptotics. Second we solve the problems for asymptotic 
coefficients. And third we determine the domain of the validity for expansion 

CD- 



3.1 The equations for coefficients 

Let us construct the internal asymptotic expansion in the domain \l\ -C 1. 
Denote 

A(xi, ti, e) = -Usxi, eti). (13) 

In the domain 1 « A « r 1 both asymptotics J2Jl and (|llfl are valid. This 
fact allows us to obtain the asymptotic representation for coefficients of the 
internal asymptotics. Substitute I = eX in formula (|10f) and expand the obtained 
expression with respect to powers of small parameter e. It yields 

n ^ — ^ . n+1 - 

Wk= 2^ x-J u {(x 2 ,t 2 ,x 1 ), ken n , x — > -00. (14) 

j=(n-fc+l) 

Let us obtain the differential equations for the coefficients of asymptotics 
l|ll|). Substitute ifTTl) . |Q in equation JIJ and collect the terms with equal 
powers of small parameter and exponents. It yields the equations for coefficients 

n 

Wk- In particularly, the terms of the order e give us the equations for the 
leading-order terms of the asympotics 

2id t2 Sd tl Wi -2id X2 Sd Xl Wi -X Wi= f, (15) 
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1 

and complex conjugated equation for W-i- 

The relation of the order e 3 in equation gives four equations. Two of 

2 2 

them are complex conjugate differential equations for Wi and W~i- 

2 2 2 1 1 

2id t2 Sd tl Wi -2id X2 Sd Xl Wi -A Wi= d 2 Xl Wi -d\ Wi - 
-i[d 2 2 S - d 2 X2 S] Wi -2id t2 Sd t2 Wi +2iO X2 Sd X2 Wi -H Wi | 2 Wi, (16) 
two another equations are algebraic. These last equations allow us to determine 

3 3 

the functions Ws and W -3 

3 7 1 , 
Ws= j(Wif. 

The higher-order terms are calculated by the same way. In particularly, 
the terms in the case lower index is equal to 1 are determined by differential 
equations. 

ti ft it it 

2id t2 Sd tl Wi ~2id X2 Sd Xl Wi -A Wi=Fi ■ (17) 
The right hand side of equation (|17J) has the form 

^2 7^ X TL 1 Th X TL X 

Fi= ~2id t2 Sd t2 W 1 +2id X2 Sd X2 W 1 + {dt 2 S) 2 W 1 -{d X2 S) 2 W 1 - 
-d 2 tl W 1 -ffig, 1 -d t2 d tl W 1 +3*^ W 1 - 

n-3 „ n— 3 « ni n 2 n 3 

-<9 t2 W 1 +9 2 2 W 1 -7 2^ W kl Wk 2 Wk 3 - (18) 

ni + n 2 + "3 = n + 1, 

fcl + &2 + &3 = 1 

kj e n nj , j = 1,2,3 

The higher-order terms in the case the lower index is not equal to 1 are deter- 
mined by algebraic equations 

n *y f n — 2 n — 2 n — 2 n — 2 

Wk= ~fc2~~[ [-Hdt 2 Sd t2 Wk+2id X2 Sd X2 Wk +(dt 2 S) 2 W k -{d X2 S) 2 W k - 

9 n — 2 n — 2 n — 3 n— 3 

-d 2 tl Wk+d 2 Xl Wk~d t2 d tl Wk+d X2 d xl Wk~ 

n-4 „ n— 4 ^-^ n t n 2 n 3 \ 

~d 2 2 W k +d 2 X2 W k - 2^ W kl Wk 2 Wk 3 J. (19) 

TT.1 + «2 + ^3 = 1 + 1) 

fci + k 2 + k 3 = k 
k 3 G 0^, j = 1,2,3 

3.2 The solvability of equations for higher-order terms 

n 

In this section we present the explicit form for higher-order term Wi and inves- 
tigate the asymptotic behaviour as A — ► ±00. 
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3.2.1 Characteristic variables 



___ 

The function Wi satisfies equation (IT7|) . The solution is constructed by char- 
acteristic method. Define the characteristic variables cr, £. We choose a point 
[x\, t\) such that d X2 l\( x o t o-j ^ as origin and denote by a the variable along the 
characteristic family for equation l|17(l . We suppose cr = on the curve A = 0. 
The variable £ mensurates the distance along the curve A = from the point 
(x1,ti). This point (x1,t1) corresponds to £ — 0. Let the positive direction 
for parameter £ coincide with positive direction of X2 in the neighborhood of 

(*?,*?)• 

The characteristic equations for l|17|l have a form 

^ = 2d t2 S(exi,eh), ^ = ~2d X2 S(ex u eh). (20) 

OCT OCT 

The initial conditions for the equations are 

a:i|a=o = a;?, ii] CT =o = «?- (21) 



Lemma 2 The Cauchy problem for characteristics has a solutions as \a\ < 
Ci£ _1 , Ci = const > 0. 

Proof. The Cauchy problem (|20|) , (|21[1 is equivalent to the system of the integral 
equations 

xi=x 1 — 2 / d X2 S(exi, eti)d(. (22) 
Jo Jo 

Substitute t 2 = (ti — t\)e, x 2 — (xi — x\)e. It yields 

pea- pea 

t 2 = 2 / d t2 S(x 2 - ex\, t 2 - et^dC, Z 2 = -2 d X2 S(x 2 - ex?, t 2 - et^dC- 
Jo Jo 

The integrands are smooth and bounded functions on the plane x 2 ,t 2 . There 
exists the constant c\ = const > such that the integral operator is contraction 
operator as e\a\ < c\. Lemma |2 is proved. 

Assumption. We assume that the change of variables — * (cr, £) is 
unique in the neighborhood of the curve A = 0. This assumption means that 
the characteristics for equation l|17|) do not touch the curve A = 0. It means 

d X2 ld X2 S - dt 2 W t2 S ^ 0. 

It is convenient to use the following asymptotic formulas for change of vari- 
ables (xi,ti) — ► (a, £). 
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Lemma 3 In the domain \a\ -C e 1 the asymptotics as e — ► o/ f/ie solutions 
for Cauchy problem $20\) . fMj) have the form 



N 



xi(a, e, e) - = -2a9 X2 5 + 2 £ eV+Wen, eti) + O^+V^ 2 ), 

N 

h(a, £, e) - t?(0 = 2a9 t2 5 + 2 ]T eV+XC^r, efc) + 0(e JV+1 tr w+2 ), 



wher 



9n = —^(9 X2 S) 



a=0 



a=0 



The lemma proves by integration by parts of equations (|22[) . 
The next proposition gives us the asymptotic formula which connects vari- 
ables a and A as a, A — > ±oo. 



Lemma 4 Let be a <C £ 1 , i/ien: 



A = 0(O<r + O(ar 2 ), 0(C) 



da 



G — > OO. 



Proof. From formula i|13|l we obtain the representation in the form 

OO 

A = X j (x 1 ,ti,e)a : '£ j ~ 1 , 

3 = 1 



where 



1 d? 

\j{xi,ti,e) = — — A(xi,ti,e)| CT=0 . 



It yields 



A = 



Let be 



dX, 

~dlr~\°-" 
dH 



a + Oiea 7 



, d 2 X . 
da 2 '' 



da 2 



> const, £ € R. 



The function dX/da is not equal to zero 

d X2 xd X2 s + d t2 xd t2 s ) = :o. 

Let us suppose dX/da > 0. It yields 



dA _ 1 

da ~ 2 



A = 0(e)<7 + O( £ a 2 ), 0(C) 



rfA 
(Zcr 



The lemma is proved. 
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3.2.2 Solutions of the equations for higher-order terms 

The higher-order terms W±i are solutions of equation (|17|) with the given 
asymptotic behaviour A — > — oo. Equation l|17|l can be written in character- 
istic variables as 

d n n n 

i-r- Wi -A Wi=Fi ■ (25) 



Lemma 5 The solution of equation |,?7| ) with the asymptotic behaviour \1J$ as 
A — > — oo has a form 



Wi= exp(-i / dCA(a;i,ti,£)) / d( Fi (xi, h, e) exp(-i / dxA(a;i, ii, e)). 

(26) 

Proof. By direct substitution we see that expression (|26|l is the solution 
of (|25|l . The asymptotics of this solution as A — > — oo can be obtained by 
integration by parts and substitution 



da 



2d t2 Sd tl -2d X2 Sd Xl . 



It yields 



n x ^ 

J=0 



( 2d t2 sd tl -2d X2 sd xl 

\ iX 



El 

iX 



X — > — oo. 



(27) 



From formula (|18|l we obtain that formulas 1)27(1 and 1(14(1 are equivalent. The 
lemma is proved. 



3.3 Asymptotics as A — > oo and domain of validity of the 
internal asymptotics 

The domain of validity of the internal expansion is determined by the asymp- 
totics of higher-order terms. In this section we show that the n— th order term 
of the asymptotic solution increases as A™ -1 when A — > oo. This increase of 
higher-order terms allows us to determine the domain of validity for internal 
asymptotics lfTT|l as A — > oo. 



3.3.1 Asymptotics of higher-order terms 

This section contains two propositions concerning asymptotic behaviour as A — > 
oo for higher-order terms in ((11(1 . The first lemma describes the asymptotic 
behaviour of higher-order terms as A — > oo and the second one contains a result 
about asymptotics of the phase function. 
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1 

Lemma 6 The asymptotic behaviour ofWi as 1 <C A <C £ has a form 

n-lj-l 



,-_n i — n \ /JO 



j=0 k=0 

{2d t2 Sd tl -2d X2 Sd x ^ j r 



£ 

0=0 



dCX(xi,tx,e)) + 

n 

Fi 



V iX 



iX 



(28) 



Proof. Let us calculate the asymptotics of the leading-order term 

l f a [<■ r 

Wi=exp(-i dC\(xi,tx,e)) d(f(xi)exp(i dxX(x 1: t 1: ej) 



exp(-£ / dC,X(xi,ti,e)) j d(f(xi)exp(i / dxX(xi,ti,e))- 
Ja J -oo Jo 

exp(-z / d(X(xi,ti,s)) / dCf{xi) exp(i / dxA(£i, ii, e))- 

JO J-CT JO 

Further by integration by parts of the last term we obtain formula l|28() as n = 1 , 
where 

W ( °' \0= daf(xi)exp(i dxX(xi,ti,e)), 

J -oo Jo 

Fi=f(xi). 

To calculate the asymptotics of Wi in formula © we use the asymptotics 

with respect to a of the leading-order term. Integral l|26[) contains the term with 

linear increase with respect to a when n = 2. We eliminate this growing part 

from integral explicitly. The residuary integral converges as a — ► oo. It can be 

l 

calculated in the same manner as it was calculated for Wi- It yields formula 
l(2"8|) as n — 2, where 

The same direct calculations are realized for the n— th order term. The lemma 
is proved. 

To complete the proof of theorem|3|we need to obtain the domain of validity 
of asymptotics Ijllf) . The formal series l|ll[) is asymptotic when 



e n+1 W 



n+l 

<Cl, e^O. 



e n W 

LemmaElgives A <C After substitution A = el we obtain I <C 1. Theorem|21 
is proved. 
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3.3.2 Asymptotics of the phase function as A — > oo 

To obtain the asymptotics as A — * oo we need to derive the asymptotics of the 
phase function in formula (|28|) . 

Lemma 7 As A — * oo: 

d£\ = ^ + - £ (d X2 S{ Xl - x?) + d t2 S(h - + 0(eX 3 ). (29) 

Proof. Substitute the asymptotics of A from lemmalHl Calculate the asymp- 
totics of the integral in formula l|29l) 



d(\(xi,ti,e) 



2 



(-d X2 ld X2 S + dt 2 W t2 S)( + 0{e?) 



2 

{-d X2 ld X2 S + 8 t2 ld t2 S)^ + 0(ea 3 ). 

The asymptotics of the phase function S(x2, t 2 ) in the neighborhood of the curve 
/ = is represented by a segment of the Taylor series. It yields 

4 = -(d X2 S( Xl - x\) + d t2 S(h - *?))+ 

e z e 

\{S X2X2 { Xl - x\) 2 + 2S X2t2 ( Xl - a;°)(ti - t°) + S t2t2 (h - t\f)+ 

o(e(l*i-*Sl + l*i-*Sl) 3 )- 

Substitute instead of {xi — x\) and (t\ — t\) their asymptotic behaviour with 
respect to e from lemma |3J This substitution and result of lemma 0] complete 
the proof of lemma [7| 

The asymptotics as A — > — 00 contains fast oscillating terms with phase 
functions kS, k G Z . The leading-order term of the asymptotics as A — > 00 
contains the oscillations with an additional phase function. We obtain this 
result from lemma El Denote this new phase function by &(x2, t2)/e 2 . The 
asymptotics of this function is obtained in lemma The nonlinearity and 
additional phase function lead to more complicated structure of the phase set 
for higher-order terms of the asymptotics as A — * 00. 

Lemma 8 The phase set K n for the n—th order term of the asymptotics as 
A — > 00 is determined by formula 

K x = ±$; K 2 = ±$,±S, K n = U il+j2+ j 3=n Xji +Xh +Xj 3 , >Xj k € K jk . 

The proof of this lemma follows from the asymptotic formula for n—th order 
term. Representation Qllfl. formula l|28|) and lemma allow us to construct the 
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asymptotics as A — + oo of the internal expansion in an explicit form 

N ,n-ln-2 
n=l ^ j=0 k=0 

exp 



i -(d X2 S(x 1 - 4) + d t2 S(h - t»)) + 0(eA d ) 



N ,00 

E £ "(E 

n=l V j=0 



2d t2 Sd tl ~2d X2 Sd Xl 
iX 



Fi 



S(h,x 2 ) 
ex P 1 4 — j — 



E e ™ ( E exp \ ik 

n=2 ^fcef2,fc^±l 



S(t 2 ,x 2 ) 



This representation and formula H19f) complete the proof of the lemma. 



(30) 



4 Post-resonance expansion 

This section contains the construction of the asymptotics of the solution for 
Q after passage through resonance. The constructed solution has the order 
e and oscillates. The envelope function of these oscillations satisfies nonlinear 
Schrodinger equation. This section consists in two parts. The first part contain 
the construction of the formal asymptotic solution. We obtain the equations for 
higher-order terms of the asymptotics. Asymptotic behaviour for higher-order 
terms as I — > —0 follows from section l3\3.2l In the second part of this section we 
determine the domain of validity for this external asymptotics near resonance 
curve l(x2, £2) = 0. The matching method gives us the initial conditions for 
higher-order terms of the asymptotics. 

The main result of this section is formulated in the following theorem. 

Theorem 4 In the domain I 3> e the formal asymptotic solution of equation 
PJ) modulo 0(e N+1 ) has a form 



N n-2 



k=0 



n . k 



U(x,t,e) = E £ "E lnfc ( £ ) ^exp{±z$(x 2 ,i 2 )/e 2 } ±*(si,ti,*a) 



±<i> 



exp{i x {x 2 ,t 2 )/e } f x {x u t 1 ,t 2 )j. (31) 



Here the function <&{x 2 ,t 2 ) satisfies eikonal equation 

(ft 2 $) 2 -(& 2 $) 2 -l = 
and initial condition on the curve 1 = 0: 



(32) 



$|;=o = S|;=o, d t2 t>i=o = d t2 S\i= . 
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The leading-order term of the asymptotics is a solution of the Cauchy problem 
for nonlinear Schrodinger equation 

1,0 „ 1,0 1,0 1,0 1,0 

2id t2 $d t2 * $ + 9| * $ + i[<9 t 2 2 $ - 94*1 * * + 7l * *l * * = 0, 

1,0 f°° f a 

® i,\l=o= daf(xi)exp(i dxX(xx,ti,e)), 
J -00 Jo 

where the £ is defined from 

= a $ ^ = a <& 

TTie coefficients v]/ ±5, are determined from Cauchy problems for linearized 

n,k 

Schrodinger equation ijffffi) . Tfte coefficients f x , ^€ fc are determined from 
algebraic equations $%yp . The set K' . = i4T ni fc\{±<i>}. 

Theorem 1 follows from theorems 2,3 and 4. 



4.1 Structure of the second external asymptotics 

Let us construct the formal asymptotic solution from theorem 0] Substitute 
1|31|) in original equation and collect the terms of the same order with respect to 
e. It yields N + 1 equations and residual of the order e N+1 . After collecting the 
terms with the same phase functions we obtain the recurrent system of equation 
for coefficients of <|31[) . 

Let us consider equations under exp(z<i>/e 2 ). The terms of the order e 1 give 
us the equation (|32[1 for the phase function of eigen oscillations. The initial data 
is determined by matching condition and represented by value of driven phase 
S on the resonance curve I = 



<H = o = s\ l=0 , d t M=o = d t2 S\ l=0 . 

The terms of the order e 2 

/ 1.0 1,0 \ 

2i ( d t2 $d tl 1< $ - d X2 $d Xl * $ I = 

give us the homogeneous transport equation 

d t2 $a tl # $ - d X2 ®d Xx i° * = 0. (33) 

This equation allows us to determine the dependence of the leading-order term 
on characteristic variable (. Equation l|33|l along the characteristics 
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can be written in the form of ordinary differential equation 

«± = 0. (35, 

1,0 

It yields vj/ $ depends on £, where the £ is defined by 

The terms of the order e 3 which oscillates as exp(z<i>/e 2 ) are 

/ 2,0 2,0 

2i I dt 2 $d tl # * - d X2 $d Xl * $ 

1,0 1,0 „ 1,0 

z[d 2 $-5 2 2 $] * $+ 7 | * *| 2 * * = 0. 

It is convenient to write this equation in the form of ordinary differential equa- 
tion in terms of characteristic variables 



2,0 



1,0 l.O^ 1 „ 1,0 

-i[d 2 t2 <$> - dl 2 <S>] * $ - 7 | * | 2 * *. (36) 



Equation (|35|> shows that the right hand side of equation l|36f) does not depend 
on £. To avoid secularities in the asymptotics we demand the right hand side of 
equation is equal to zero. It allows to determine the dependence of the leading- 
order term on slow variable i 2 

2id t2 <S>d t2 ¥ $ + p tl 2 - (d Xl 2 ]dh * * + 

+i[dl<P - d 2 X2 <S>] * $ + 7 | * $| 2 * $ = 0. (37) 
The equations for the higher-order terms are obtained by the same manner 

(n+l,fc n+l.k \ n,fc n,/c 

<9t 2 $d tl * $-d X2 $d Xl * $1 =2idt 2 $d t2 * * - * *- 

-i[d 2 2 $-<9 2 2 $] * $+a tl ^| t2 * $-7 ^ * Q * p * 4l 

ki,k 2 ,li,l 2 ,mi ,m2 

where fei + Zi + mi = n + 2, fe 2 + Z2 + m 2 — k, a + (3 + 6 — <f>, a G Kk lt k 2 > & & 
Ki t .i 2 , 8 S K mi}m2 . 
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To construct the uniform asymptotic expansion with respect to ( we obtain 
the linearized Schrodinger equation for higher-order term 

71, fe n n,k n,k 

2id t2 $<9 t2 * 4> + <9f € * * -H[<9 t 2 2 <J>-<^ 2 $] * $ = 

n— l,fc x — ^ fci,fc2 ill, 22 mi, ma 

-d tl £dl t2 f*- 7 2^ (38) 

ki,k 2 ,li,l 2 ,mi ,fi2 ,ck,/3,<5 

where fci + ii + mi = n + 2, k 2 + h + m.2 = k, a + (3 + S = <f>, a G Kk lt k 2 , P G 
Ki u i 2 , 5 € i^ mi: m 2 . 

n 

The amplitudes \I> X as x ^ ±3> arc determined by algebraic equations 

[-(Xt 2 ) 2 + (Xx 2 ) 2 + l] *x=F x , X^±1». (39) 

Here the right hand side of the equation depends on previous terms and their 
derivatives 

n,k n—l.k n—l.k n— 2,fc n—2.k 

F x =-2i X t 2 d tl * x +2ixx 2 d xi * x - 2i Xt 2 d t2 * x -i [xt 2 t 2 - Xx 2 x 2 ] * x - 

n— 3,fc n— 4,fc ^ — T feijfea 2i^2 mi,m2 

5 txt 2 * X ~^ a fa *X-7 2^ * a * /3 # ,5, 

fci,fc2,^i^2,"H ,m2 ,a,/3,<5 

where fei + Zi + mi = n — 4, k 2 + l 2 + m 2 = k, a + (3 + 8 = x, ot G K ku k 2 , & G 

K tl j 2 , 6 G ifmi,ma< 

These equations are similar to equations for amplitudes from pre-resonance 
section. 

The obtained result is formulated below 

Lemma 9 The coefficients of formal asymptotic solution QUlp satisfy recurrent 
system of equations f3gp. (gffi , Ell- 

The right hand side of equation (|38|l has a singularity as I — > 0. The singu- 

larity appears due to vj/ x as x The analysis of the right hand side of 

the equation allows us to calculate the order of singularity as I — » 0. It is equal 
to 0(l-( n ~^). Below we prove the solvability of equation l|38|l with the given 
asymptotics as I — > 0. 

Lemma 10 The asymptotics as I — > o/ £/ie solution of equation fffffi) has the 
form 

k 1 j ~ l k 

* *(«i,ti,*2) = E E * • Tn («i»*i) ^(lnir + OW. 1-0.(40) 

j=-(n-2) "i=0 

Proof. Determine the order of the singularity of the right hand side of the 
equation as I — > 0. First consider equation l|38|) for n = 3, fc = 0. The solution 

3,0 

of this equation gives us the coefficient if $ . The nonlinearity contains the term 
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2,0 1,0 2,0 

I s 1 2 ^<s>- The function vj/g has the singularity of the order l^ 1 as I — > 0. 
It determines the order of singularity for right hand side l~ 2 . We construct the 

3,0 

asymptotics of 4> x m ^ ne f° rm 

3,0 3,0 , n , 3,0 n . 3,0,, 3,0 

* *=* I ' * «" M0+ * » nn(0+ (41) 

Substitute jlT} in equation for n = 3. It leads to recurrent system of equations 

3 '° a k) 

for coefficients * $ 

3,0 , , m 1,0 2,0 „ „ 

-2id t2 <t>d t2 i * £ m = - I * s\ 2 i 2 , 

2id t2 $d t2 l * £' X) = L[4 <j°' X) ]. 
Here we denote the linear operator by 



3,0 

The regular part \j> $ of the asymptotics satisfies the nonhomogeneous linear 
Schrodinger equation. The right hand side of the equation is smooth 

3,0 3,0 n 3,0 f , u 

L[* 9 ] = -lhi.\l\L[* ^ 1) ]-2id t2 ^d t2 l * £• ■ 

The initial condition for the regular part of the asymptotics is determined below 
by matching with the internal asymptotic expansion. 

n.k 

The structure of the terms ±$ for n > 3 has a similar form. The right 
hand side of equation (|38|l depends on junior terms. These singularities can be 
eliminate 

-(n-2) -j+l , —7 

n.k i — , . . n,k , . s n.k 

j—0 m—0 

n,k , . . 

The coefficients / ^ do not contain singularities as I — > 0. These coefficients 
are easy calculated. 

The direct substitution of (|4t)|l in equation and collecting the terms with the 
same order of I complete the proof of lemma HHI 

4.2 The domain of validity of the second external asymp- 
totics and matching procedure 

The domain of validity of the second external asymptotics is determined by 

n+l 

V 
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Formulas (|31[l and l|4(J|) give the condition 

I < e. 

The domain \l\ <C 1 of validity of the internal asymptotics and domain of 
validity of the second external asymptotics are intersected. This fact allows 
to complete the construction of the second external asymptotics by matching 
method |17| . The structure of singular parts of the internal asymptotics as 
A — > +00 and external asymptotics as / — > are equivalent. The coefficients 
are coincided due to our constructions. The matching of regular parts of these 
asymptotics takes place due to 

\ l=0 =W {m (0- 

The function W (°>°)(£) is determined in lemma El 

In particular, the initial condition for the leading-order term has a form 

1,0 f°° f 7 

<I> $|;=o = / daf{xi) exp(i / dxX(xi,t 1 ,e)). 
J -00 Jo 

The soliton theory for nonlinear Schrodingcr equation leads us to the fact that 

1,0 

the function vj/ $ contains the solitary waves when f(x\) is sufficiently large. 
Theorem 0] is proved. 

Acknowledgments. We are grateful to I.V. Barashenkov, L.A. Kalyakin 
and B.I.Sulcimanov for helpful comments and for help in improving of the math- 
ematical presentation the results. 
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